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Pressure  surge  analysis  of  petroleum  pipeline  transporting 
automotive  gas  oil  (AGO)  also  known  as  Diesel  oil  was  carried  out  in 
this  research  work.  Pressure  transient  analysis  is  often  more 
significant  than  the  steady  state  analysis  that  engineers  usually  use  in 
pipeline  design.  Pressure  transient  analysis  helps  to  understand  the 
additional  pressures  the  pipeline  can  be  subjected  to  as  a  result  of 
instantaneous  rapid  valve  closures  or  pump  failure.  The  fluid  pressure 
and  flow  rate  in  the  pipeline  system  may  change  significantly  at  some 
intervals  of  time  due  to  the  valve  closure  and  such  types  of  unsteady 
situations  are  encountered  more  often  in  pipelines  where  the  valves 
are  suddenly  closed.  In  this  paper,  pressure  surge  due  to 
instantaneous  valve  closure  in  a  petroleum  pipeline  conveying  AGO 
was  studied  in  a  virtual  environment.  WANDA  Transient  4.5.1210 
commercial  software  was  used  for  the  analysis  of  the  pressure  surge 
in  the  pipeline  due  to  instantaneous  valve  closure  time  of  4.75s.  It  was 
observed  in  the  study  that  pressure  at  some  nodes  rise  significantly  up 
to  about  1400  kPa  against  the  initial  inlet  pressure  of  120  kPa  due  to 
the  instantaneous  valve  closure  and  it  drastically  drops  at  some  nodes 
to  negative  pressure  of  about  -100  kPa  and  hence  the  formation  of 
cavitations.  The  analysis  showed  that  the  magnitude  of  the  pressure 
surge  decreases  as  the  valve  closure  is  increased. 


1. Introduction 

Pressure  surge  also  termed  pressure  transient,  hydraulic  transient,  transient  waves,  fast  transients, 
fluid  transients,  hydraulic  hammer,  oil  hammer  or  water  hammer  refers  to  rapid  changes  of  pressure 
due  to  changes  in  some  of  the  flow  parameters  in  a  pipe  system  [1,2].  Pressure  surge  has  destructive 
and  catastrophic  consequences,  such  as  collapsing  pipes  and  ruptured  valves  [1,  3].  [4]  Reported 
that  there  are  two  categories  of  damage  caused  by  pressure  transient  events.  Namely;  catastrophic 
failure  and  fatigue  like  failure.  Catastrophic  failure  is  a  type  of  failure  caused  by  high  magnitude 
transient  waves  generated  as  a  result  of  valve  closure  or  pump  failure  while  a  fatigue-like  failures 
are  normally  caused  as  a  result  of  prolong  repeated  impacts  of  smaller  magnitude  transient  pressure 
over  a  long  period.  Therefore,  it  is  imperative  to  understand  and  investigate  the  phenomena  and  its 
causes  as  well  as  accurately  calculate  and  analyze  its  effects  on  the  pipeline.  The  importance  of 
carrying  out  pressure  surge  analysis  cannot  be  over  emphasised,  transients  analyses  are  carried  out 
in  pipelines  networks  in  order  to  verify  whether  the  networks  are  operating  within  acceptable 
maximum  operating  pressure  as  well  meets  the  regulations  and  standards  [5]. 
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[6]  Reported  that  pressure  surge  analysis  and  evaluation  play  a  very  vital  role  in  the  design,  operation 
and  maintenance  of  new  and  existing  pressurized  pipeline  systems.  Also  [6  and  7]  stated  that 
pressure  surge  occurs  in  a  pipeline  network  as  a  result  of  either  pump  failure  or  due  to  sudden  valves 
closure  or  opening.  This  paper  is  aimed  at  investigating  the  pressure  surge  due  to  instantaneous 
valve  closure  in  a  petroleum  pipeline  conveying  automotive  gas  oil  (AGO)  using  simulation 
approach.  Conducting  transient  analysis  in  a  pipeline  system  is  often  more  important  than 
conducting  steady  state  condition  analysis  in  a  pipeline  [6].  There  are  quite  large  number  researches 
conducted  on  hydraulic  transients  occurring  in  closed  conduits  [8,  9,  10,  11,  12,  and  13].  Pressure 
surge  analysis  methods  range  from  analytical  methods  to  numerical  solutions  [9].  According  to  [14] 
these  methods  are  further  divided  into  either  elastic  or  rigid  column  method.  Elastic  method  is  a 
method  of  transient  analysis  that  involves  solving  partial  differential  equations.  Elastic  method  also 
involves  evaluating  the  acoustic  pressure  wave.  While  a  rigid  column  method  is  a  method  of 
pressure  surge  analysis  that  involves  solving  simple  ordinary  differential  equations  mathematically 
or  numerically.  In  this  method,  the  elasticity  of  the  pipe  and  the  compressibility  of  the  fluid  are 
ignored  in  the  analysis  and  whole  of  the  fluid’s  column  is  assumed  to  move  as  a  rigid  body  [15].  In 
both  cases,  quasi-linear  hyperbolic  partial  differential  equations  are  used  in  the  analysis  of  unstable 
fluid  flow  in  pipelines  [16]. 

Some  of  the  methods  used  in  pressure  surge  analysis  are  arithmetic  mean  method  [17],  Graphical 
method  [18],  analytical  [19],  experimental  [20],  Method  of  characteristics  (MOC)  [21],  Einite 
difference  methods  (EDM)  [22],  Wave  plan  method  [23  and  24].  [10]  reported  that  the  most  widely 
accepted  and  used  methods  of  pressure  surge  analysis  are  the  method  of  characteristics  (MOC)  and 
wave  characteristics  method  (WCM)  and  the  main  distinction  between  the  two  methods  is  the  way 
pressure  waves  are  traced  between  pipe  boundaries.  The  MOC  use  numerical  method  to  trace  a 
disturbance  in  a  grid  on  characteristics,  whereas  WCM  uses  wave  propagation  method  to  trace  the 
disturbance.  These  two  outstanding  methods  are  well  documented  in  pressure  transient  [16  and  25] 
and  have  been  implemented  in  various  computer  programs  for  pipeline  system  transient  analysis. 
Nowadays  computational  fluid  dynamics  (CED)  method  is  used  to  carry  out  fluid  flow  and  flow 
condition  analysis.  In  this  research  work,  CED  simulation  software  called  WANDA  Transient 
4.5.1210  was  used  for  the  investigation  of  pressure  surge  in  a  pipeline  network.  WANDA  is  one  of 
the  most  outstanding  commercial  simulation  software  that  uses  MOC  for  the  analysis  of  fluid  and 
heat  flow  in  pipeline  networks  [26] . 


2.  Methodology 

2.1  Data  Collection 

The  fluids  and  pipeline  parameters  required  for  the  analysis  as  presented  in  Tables  1  and  2  were 
obtained  from  direct  measurement,  from  the  archives  of  the  Nigeria  National  Petroleum  Company 
Eimited  (NNPC)  and  literature. 


Table  1:  Pipe  parameters 


Pipe 

Pipe 

Diameter 

Thickness 

Roughness  (e/d) 

Ey 

Material 

(m) 

(m) 

(mm) 

(N/mh 

1 

Carbon  steel 

0.3556 

0.016 

0.045 

210x10'^ 
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Table  2:  Fluid  parameters 

Fluid  Density  Flow  rate  Pressure 

(mVkg)  Q  (m^/s)  (kN/m^) 

AGO  870  0.297  120 

Change  in  pressure  of  fluid  flowing  in  a  closed 
change.  The  basic  water  hammer  equation  is  used 
by  the  surge  in  the  pipes,  as  shown  in  Equation  1. 


F 

Ey 

Cp 

K 

(N/m2) 

(N/m2) 

(kJ/kg/K) 

(N/m2) 

0.006 

213.84 

2.22 

1.07x10® 

conduit 

is  directly  proportional  to  the  velocity 

to  express  the  change  in  pressure  head  produced 


AH  =  --Av  (1) 

a 

[21]  also  reported  that  the  fundamental  water  hammer  theory  which  describes  the  pressure  amplitude 
was  laid  by  Joukowsky  and  is  presented  in  Equation  2. 


AP  =  ±paAV  or  AH ±=^  (2) 

Equation  (3)  was  used  for  calculating  wave  speed  propagation  for  transient  flow  in  the  pipeline  as 
reported  by  [6  and  27] . 


a  = 


k 


(3) 


According  to  [7]  the  bulk  modulus  of  elasticity  of  a  fluid  is  an  important  parameter  in  the  analysis 
of  wave  speed  of  fluid  and  it  can  be  obtained  by  using  Equation  4. 


K  = 


AP 

Ap/p 


(4) 


The  head  loss  due  to  water  hammer  in  the  pipeline  can  be  calculated  using  the  Darcy-Welsbach 
relation  presented  as  Equation  5: 


hf  = 


fDLV^ 

2gD 


(5) 


The  relative  roughness  of  the  pipe  is  calculated  using  Equation  6, 


Relative  roughness  = 


(6) 


[28]  Reported  that  instantaneous  valve  closure  is  characterized  by  valve  closure  time  less  than  T 
and  the  value  of  T  can  be  calculated  by  using  Equation  7. 
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T 


2L 


a 


(7) 


HT  behaviours  in  closed  conduits  can  be  analyzed  by  using  Equations  of  motion  and  continuity  [29], 
the  Equation  is  shown  in  Equation  8  as  reported  by  [27] 


dV 


4-  T/ 


dV  ,  1  dp 


(8) 


The  term  VdV/dx  in  Equation  8  is  neglected  in  transient  analysis  as  a  result  of  low  Mach-number 
and  unsteady  flows.  [7]  Reported  that  Equation  8  will  reduce  to  Equation  9 


dV  1  dp  ^  ,  fV\V\ 

- h  -f  g  Sin  0  -I - = 


n 


rii- 


(9) 


Also  according  to  [29]  the  general  form  of  the  continuity  Equation  can  be  presented  in  Equation  10. 


1  dp  ^  1  dA  ^  dV 
p  dt  A  dt  dx 


(10) 


If  an  elastic  pipe  is  filled  with  a  compressible  fluid,  Equation  10  will  reduce  to  relation  presented  in 
Equation  11. 


dp  dv  2^^  n 

—  +  V - h  pa^  —  =  0 

dt  dx  ox 


(11) 


The  friction  factorcan  be  calculated  by  using  the  Colebrook- White  Equation  as  shown  in  Equation 
12  [30] 


// 


-21ogio 


H  2.51 
3.7D  Re^f 


2.2  Mechanics  of  Water  Hammering 


(12) 


In  this  study,  a  hypothetical  petroleum  pipeline  network  is  adopted  for  this  research  work.  Data  used 
in  this  research  were  obtained  from  literature,  NNPC  archives  and  field.  The  pipeline  network 
consists  of  the  followings,  upstream  and  downstream  reservoirs,  pipes,  pump  and  valves  (gate  valve 
and  non-return  valves).  The  pipeline  networks  parameters  under  study  are:  upstream  reservoir  (Bl) 
at  a  surface  elevation  of  14.1m  with  a  pressure  of  120  kPa.  The  pipeline  network  is  made  up  of  four 
carbon  steel  pipes  of  equal  diameters  of  0.3556m,  thickness  of  0.016m  and  surface  roughness  of 
0.045mm.  All  the  pipes  are  connected  in  series  between  the  upstream  reservoir  (Bl)  and  the 
downstream  reservoir  (B2).  The  upstream  end  of  the  first  pipe.  Pipe  1  (PI)  is  connected  to  the 
upstream  reservoir  via  node  A  while  its  downstream  end  is  connected  to  the  pump  at  node  B,  pipe 
PI  has  a  length  of  1000m.  Pipe  2  (P2)  is  the  second  pipe  in  the  series,  its  upstream  end  is  connected 
to  the  pump  at  node  C  and  its  downstream  is  connected  to  the  check  valve  at  node  D,  it  has  a  length 
of  10000m.  The  third  pipe  in  the  series  of  the  pipes  is  pipe  3  (P3),  P3  is  also  connected  the  check 
valve  at  its  upstream  at  node  E  while  at  the  downstream  is  connected  to  the  gate  valve  at  node  E;  P3 
also  has  a  length  of  10000  m.  The  last  pipe  in  the  series  is  pipe  4  (P),  it  has  a  length  of  1000m,  the 
upstream  end  of  P4  is  connected  to  the  gate  valve  at  node  G  while  its  downstream  end  is  connected 
to  the  downstream  reservoir  of  the  pipeline  network  at  node  H.  The  valve  closure  times  used  in  this 
analysis  are  4.75  s,  9.5  s  and  19  s  respectively. 
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B  C  D  E  F  G 

P1  P1  P2  C1  P3  F1  P4 


Figure  1:  Pipeline  Network  System 

In  a  fluid  transporting  pipeline  network,  if  a  valve  is  closed  instantaneously,  the  momentum  of  the 
fluid  will  be  shattered  and  high  pressure  wave  will  be  built-up  accordingly.  This  built-up  high 
pressure  wave  will  be  transmitted  down  the  pipe  length  with  the  velocity  of  the  sound  wave  that 
may  leads  to  knocking. 

Table  3  presents  some  of  the  parameters  inputted  into  the  WANDA  software  for  the  analysis  of  the 
hydraulic  transient  in  AGO  while  Table  4  presents  some  physical  constants  under  which  the 
simulation  was  conducted. 


Table  3:  Properties  of  AGO 
AGO 


Liquid  name 
Rheology  type 
Density 
Bulk  modulus 
Vapour  pressure 
Kinematic  viscosity 


Newtonian 
871.0  kg/m^ 

1. 477x1  O'*  N/m2 
1.100  kPa 
2.860e-5  mVs 


Table  4:  Physical  constant  used  in  the  simulations 


Atmospheric  pressure 

101.4  (kPa) 

Gravitational  acceleration 

9.810  (m/s2) 

Ambient  temperature 

37.00  (°C) 

3.  Results  and  Discussion 

Figure  2  shows  the  pressure  head  distribution  before  the  advent  of  the  valve  closure  in  the  pipeline 
network.  That  figure  indicated  that  pressure  head  starts  to  rise  from  node  B  and  reached  a  peak  of 
0. 107164  m  at  node  C,  this  is  due  to  the  presence  of  pump  at  that  node  which  adds  more  pressure  to 
the  fluid.  After  node  C,  the  pressure  starts  to  drop  across  nodes  D,  E  and  F.  The  pressure  stabilises 
between  nodes  G  and  H  a  pressure  head  of  0.3  m.  It  was  observed  in  this  study  that  the  pressure 
head  was  maximum  at  the  upstream  of  the  pump  after  development  of  first  pressure  wave  which  is 
due  to  instantaneous  closure  of  the  valve  at  the  downstream  end  of  the  pipe.  This  condition  is  in 
agreement  with  what  reported  by  [6].  The  pressure  head  was  not  dampened  but  they  were  found  to 
be  oscillating  by  increasing  and  decreasing  in  the  total  simulation  time  which  may  be  due  to  wave 
propagation  time  and  consideration  of  series  pipe  in  the  current  study. 
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0.12 
^  0.1 
I  0.08 


ABCDEFGH 


Hydraulic  Nodes 


■  AGO 


Figure  2:  Pressure  heads  at  various  hydraulic  nodes  of  the  pipeline 

3.1  Simulated  results  for  AGO  at  instantaneous  valve  closure  times  of  4.75  s,  9.5  s  and  19  s. 

The  effects  of  the  instantaneous  valve  closure  were  observed  as  pressure  oscillation  at  different 
types  of  closing  times  as  reported  by  [32].  The  results  of  the  analysis  are  in  agreement  with  [33] 
which  reports  that  the  pressure  wave  will  travel  back  and  forth  in  the  pipeline  until  the  kinetic  energy 
is  dissipated  by  friction.  And  this  process  will  occur  both  upstream  and  downstream  from  the  valve. 
However,  the  initial  pressure  will  increase  on  the  upstream  side  of  the  valve  and  decrease  on  the 
downstream  of  the  valve  as  shown  in  figures  3,  4  and  5. 

3.1.1  Simulated  results  for  AGO  at  instantaneous  valve  closure  times  of  4.75  s 

Figure  3  present  results  of  pressure  transients  due  to  instantaneous  valve  closure  time  of  4.75s.  The 
figure  depicts  that  there  is  variation  in  pressure  in  the  various  hydraulic  nodes  of  the  pipeline.  The 
pressure  of  the  fluid  under  normal  operating  conditions  at  nodes  A,  B,  C,  D,  E,  F,  G  and  H  are 
118.72  kPa,  60  kPa,  1000  kPa,  540  kPa,  560  kPa,  470  kPa,  260  kPa  and  254.73  kPa  respectively. 
But  when  the  gate  valve  is  closed  instantaneously,  the  pressure  of  the  fluid  fluctuates  and  oscillates 
between  minimum  and  maximum  pressures  before  stabilising  as  shown  in  Figure  3. 

The  pressure  at  node  A  fluctuates  and  oscillates  between  minimum  pressure  of  1 18.4  kPa  and  1 18.9 
kPa  at  times  of  5  s  and  68.1  s  respectively.  At  node  B,  the  pressure  of  the  fluid  oscillates  between 
-  100.5  kPa  and  367  kPa.  The  development  of  negative  pressure  at  node  B  also  leads  to  the  formation 
of  cavitation  voids  at  the  node  and  in  pipe  PI  respectively.  At  node  C  the  pressure  of  the  fluid 
fluctuates  between  347.4  kPa  and  1262  kPa.  The  pressure  of  the  fluid  drops  to  a  minimum  pressure 
of  295.5  kPa  at  a  time  of  8.3  s  before  rising  to  a  maximum  pressure  of  1307  kPa  after  59.1  s.  The 
Figure  also  depicts  that  pressures  at  node  E  and  F  fluctuates  and  oscillates  between  minimum 
pressures  of  295.5  kPa,  223.7  kPa  at  times  of  9.8  s  and  16.5  s  and  maximum  pressures  of  1317  kPa 
and  1417  kPa  respectively.  Nodes  G  and  H  are  the  two  nodes  at  the  downstream  side  of  the  gate 
valve.  The  pressure  of  the  fluid  at  the  nodes  also  fluctuates  between  minimum  pressures  of  106.8 
kPa  at  16.5  s  and  254.8  kPa  occurring  at  187.4  s  respectively. 
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AGO  4 


Tme  (s) 

—  Pressure  H-node  A  —  Pressure  H-no<Je  B  —  Pressure  H-node  C  —  Pressure  H-node  D 

—  Pressure  H-node  E  —  Pressure  H-node  F  —  Pressure  H-node  G  —  Pressure  H-node  H 

Figure  3:  Pressure  Transients  at  Hydraulic  Nodes  for  a  valve  closure  time  of  4.75  s 

3.1.2  Simulated  results  for  AGO  at  instantaneous  valve  closure  times  of  9.5  s 

Figure  4  present  results  of  pressure  transients  due  to  instantaneous  valve  closure  time  of  9.5  s  in  the 
pipeline  network.  The  figure  depicts  that  there  is  variation  in  pressure  in  the  various  hydraulic  nodes 
of  the  pipeline.  The  pressure  of  the  fluid  under  normal  operating  conditions  at  nodes  A,  B,  C,  D,  E, 
F,  G  and  H  are  1 18.72  kPa,  60  kPa,  1000  kPa,  540  kPa,  560  kPa,  470  kPa,  260  kPa  and  254.73  kPa 
respectively.  But  when  the  gate  valve  is  closed  instantaneously,  the  pressure  of  the  fluid  fluctuates 
and  oscillates  between  minimum  and  maximum  pressures  before  stabilising  as  shown  in  Figure  4. 

The  pressure  of  the  fluid  fluctuates  and  oscillates  between  minimum  and  maximum  pressures  as  a 
result  of  the  instantaneous  valve  closure.  The  pressure  at  node  A  fluctuates  and  oscillates  between 
minimum  pressure  of  1 18.4  kPa  and  1 19  kPa  at  times  of  5  s  and  28.5  s  respectively.  At  node  B,  the 
pressure  of  the  fluid  rises  to  376  kPa  at  a  time  of  0.1  s  but  later  after  about  3.5  s  the  pressure  drops 
to  -100.5  kPa.  The  pressure  keeps  oscillating  before  stabilisation.  The  drop  in  pressure  to  a  negative 
pressure  leads  to  the  development  of  cavitation  voids  at  the  node  and  along  pipe  PI.  At  node  C  the 
pressure  of  the  fluid  fluctuates  between  a  minimum  pressure  value  of  347.4  kPa  and  a  maximum 
pressure  of  1264  kPa  at  times  of  0.1  s  and  70.5  s  respectively.  The  pressure  of  the  fluid  at  nodes  D 
and  E  drops  to  a  minimum  pressure  of  295.5  kPa  at  a  time  of  8.3  s  and  9.3  s  respectively  but  later 
the  pressure  at  D  rise  to  a  maximum  pressure  of  1303  at  a  time  of  59.1  s.  Also,  the  pressure  at  node 
E  rises  to  1307  kPa  at  a  time  of  131.2  s.  The  figure  also  depicts  that  pressure  at  node  F  fluctuates 
and  oscillates  between  minimum  pressures  of  200. 1  kPa  at  a  time  of  16.5  s  and  maximum  pressures 
of  1391  kPa  at  a  time  of  18. 1  s.  Nodes  G  and  H  are  the  two  nodes  at  the  downstream  side  of  the  gate 
valve.  The  pressure  of  the  fluid  at  the  nodes  fluctuates  between  minimum  pressures  of  101.8  kPa  at 
16.5  s  and  254.6  kPa  occurring  at  187.4  s  and  a  maximum  pressure  of  565.8  kPa  and  254.8  kPa 
respectively  at  times  of  18.9  and  197.4s. 
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—  Pressure  H-no<3e  A  —  Pressure  H-no<Je  B  —  Pressure  H-no<Je  C  —  Pressure  H-ncxJe  D 

—  Pressure  H-node  E  —  Pressure  H-no<Je  F  —  Pressure  H-node  G  —  Pressure  H-nude  H 


Figure  4:  Pressure  Transients  at  Hydraulic  Nodes  for  a  valve  closure  time  of  9.5  s 

3.1.3  Simulated  results  for  AGO  at  instantaneous  valve  closure  times  of  19  s 

Figure  5  present  results  of  pressure  transients  due  to  instantaneous  valve  closure  time  of  19  s  in  the 
pipeline  network.  The  Figure  depicts  that  there  is  variation  in  pressure  in  the  various  hydraulic  nodes 
of  the  pipeline.  The  pressure  of  the  fluid  under  normal  operating  conditions  at  nodes  A,  B,  C,  D,  E, 
F,  G  and  H  are  1 18.72  kPa,  60  kPa,  1000  kPa,  540  kPa,  560  kPa,  470  kPa,  260  kPa  and  254.73  kPa 
respectively.  But  when  the  gate  valve  is  closed  instantaneously,  the  pressure  of  the  fluid  fluctuates 
and  oscillates  between  minimum  and  maximum  pressures  before  stabilising  as  shown  in  Figure  5. 

The  pressure  of  the  fluid  fluctuates  and  oscillates  between  minimum  and  maximum  pressures  as  a 
result  of  the  instantaneous  valve  closure.  The  pressure  of  the  fluid  at  node  A  fluctuates  and  oscillates 
between  minimum  pressure  of  1 18.4  kPa  and  1 19  kPa  at  times  of  5  s  and  134.2  s  respectively.  At 
node  B,  the  pressure  of  the  fluid  rises  to  367  kPa  at  a  time  of  0.1  s  but  later  after  about  3.5  s  the 
pressure  drops  to  -100.5  kPa.  The  drop  in  pressure  to  a  negative  pressure  leads  to  the  development 
of  cavitation  voids  at  the  node  and  along  pipe  PI.  At  node  C  the  pressure  of  the  fluid  fluctuates 
between  a  minimum  pressure  value  of  347.4  kPa  and  a  maximum  pressure  of  1260  kPa  at  times  of 
0.1  s  and  139.3  s  respectively.  The  pressure  of  the  fluid  at  nodes  D  and  E  drops  to  a  minimum 
pressure  of  295.5  kPa  at  a  time  of  8.3  s  but  later  the  pressures  at  these  two  nodes  rise  to  a  maximum 
pressure  of  1293  kPa  at  a  time  of  134.9  s.  The  figure  also  depicts  that  pressure  at  node  E  fluctuates 
and  oscillates  between  minimum  pressures  of  180.6  kPa  at  a  time  of  16.5  s  and  maximum  pressures 
of  1371  kPa  at  a  time  of  57.9  s. 

Nodes  G  and  H  are  the  two  nodes  at  the  downstream  side  of  the  gate  valve.  The  pressure  of  the  fluid 
at  nodes  G  fluctuates  between  minimum  pressures  of  94.94  kPa  at  16.5  s  and  maximum  pressure  of 
578.6  kPa  occurring  at  18. 1  s  while  the  pressures  at  node  H  fluctuates  between  a  maximum  pressure 
of  254.8  kPa  and  254.4  kPa  respectively  at  times  of  18.9  and  161.1s. 
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—  Pressure  H-no<Je  A  —  Pressure  H-node  B  —  Pressure  H-node  C  — -  Pressure  H-node  D 
^  Pressure  H-node  E  ^  Pressure  H-node  F  ^  Pressure  H-node  G  ^  Pressure  H-node  H 

Figure  5:  Pressure  Transients  at  Hydraulic  Nodes  for  a  valve  closure  time  of  19  s 

3.2  Comparison  of  the  high-pressure  transients  due  to  various  valve  closure  times 

The  comparison  of  the  simulated  results  for  various  valve  closure  times  in  the  pipeline  transporting 
the  AGO  is  shown  in  Figures  6.  The  graph  depicts  that  pressure  transient  due  to  instantaneous 
closure  of  a  gate  valve  in  a  petroleum  pipeline  transporting  AGO  reduces  as  the  valve  closure  time 
increases.  Therefore,  in  a  pipeline  network,  it  is  better  to  have  a  longer  valve  closure  time  so  as  to 
reduce  or  eliminate  the  possibility  of  pressure  transients,  column  separation  or  the  formation  of 
cavitation  voids  that  may  eventually  leads  to  pipeline  failure.  The  results  obtained  in  this  analysis 
are  in  agreement  with  the  results  reported  by  [31]  in  which  the  rise  in  pressure  inside  network  has 
different  magnitudes  depending  upon  the  valve  closure  times.  The  valve  closure  time  has  to  be 
sufficiently  increased  to  avoid  drastic  pressure  rise  inside  the  pipe  due  to  water  hammer 
phenomenon. 
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Figure  6:  Comparison  of  pressure  transient  in  AGO  at  different  valve  closure  times 
4.  Conclusion 

Effects  of  hydraulic  transient  due  to  instantaneous  valve  closure  in  a  petroleum  pipeline  were 
investigated  in  this  study.  It  was  observed  in  the  study  that  wave  and  pressure  propagations  in  the 
pipeline  network  oscillate  between  high-  and  low-pressure  values  due  to  valve  closure.  The  study 


98 


A.  B.  Muhammad  et  al./  NIPES  Journal  of  Science  and  Technology  Research 

2(2)  2020  pp.  90-101 

showed  that  the  velocities  at  nodes  after  the  gate  valve  reduce  to  zero.  Pressure  rise  was  also 
observed  at  the  node  where  the  valves  are  located.  In  this  research  WANDA  Transient  4.5.1210 
simulation  software  was  used  to  analyze  pressure  surge  in  the  pipeline  network  due  the 
instantaneous  closure.  Pressure  fluctuations  were  observed  in  the  simulation  analysis  as  a  result  of 
the  valve  closure  and  nodes  F  where  the  gate  valve  is  located  records  highest  pressure  surge  while 
nodes  B  a  node  before  the  pump  records  negative  pressures  and  cavitations  in  all  the  pipeline 
networks  considered  with  different  valve  closure  time.  The  research  found  that  as  the  time  of  valve 
closure  increases,  the  magnitude  of  pressure  surge  developed  decreases.  The  research  recommends 
that  surge  tank  should  be  installed  at  node  F  to  stabilize  the  pressure  surge  and  also  air  vessels  are 
to  be  installed  at  nodes  B  to  curtail  damages  due  to  cavitations  at  the  node  and  pipe  PI. 

Notations 

A  Area  of  pipe,  m? 

C  Shock  wave  speed,  m/s 

Cp  Specific  heat  of  fluid  at  constant  pressure,  kJ/kgK 

d  Diameter  of  pipe,  m 

Ey  Young’s  Modulus  for  the  pipe  material,  N/m^ 

f  Moody  friction  factor 

g  Acceleration  due  to  gravity,  m/s^ 

hf  head  loss  due  to  friction,  m 

K  Fluid  bulk  Modulus,  N/m^ 

L  Length  of  pipe,  m 

m  Mass  flow  rate.  Kg/s 

P  Pipe  pressure,  N/m? 

Po  Stagnation  pressure,  N/m? 

Re  Reynolds  number 

R  resulting  force  on  the  bend  (N) 

Rx  resulting  force  in  x-direction  (N) 

Ry  resulting  force  in  x-direction  (N 

t  Period  (time),  s 

th  Pipe  thickness,  m 

T  Temperature,  K 

To  Stagnation  temperature,  K 

V  Flow  velocity,  m/s 

Greek  Letters 

p  Mass  density  of  fluid,  Kg/m^ 

y  Ratio  of  specific  heats 

p.  Dynamic  viscosity  (N/m^) 

71  pye 
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